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Abstract 

Recently Gambini and Pullin proposed a new consistent discrete approach to quantum gravity 
and applied it to cosmological models. One remarkable result of this approach is that the cosmo- 
logical singularity can be avoided in a general fashion. However, whether the continuum limit of 
such discretized theories exists is model dependent. In the case of massless scalar field coupled 
to gravity with A = 0, the continuum limit can only be achieved by fine tuning the recurrence 
constant. We regard this failure as the implication that cosmological constant should vary with 
time. For this reason we replace the massless scalar field by Chaplygin gas which may contribute 
an effective cosmological constant term with the evolution of the universe. It turns out that the 
continuum limit can be reached in this case indeed. 
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I. INTRODUCTION 



The cosmological singularity is a longstanding and fundamental problem in quantum 
cosmology. Nowadays both string theory and loop quantum gravity have been taking this 
issue into account seriously, see [J Q and for recent progress along both approaches 
respectively. Motivated from string theory, the cosmological solution from big crunch to big 
bang has been proposed in [lj, which can also be viewed as the orbifold compactification 
of solutions in higher dimensional field theory. This solution has very peculiar features. 
For instance it can be viewed as the Milne universe but identified by the boost such that 
it contains closed timelike loops and more badly, not Hausdorff. On the other hand this 
solution plays a special role in string/M theory since such spacetime is flat everywhere 
except at the singularity and the scalar field keeps staying in the weak coupling region while 
approaching to the singularity. Consequently the behavior of strings nearby such singularity 
in context of string/M theory has been studied. Although some positive results have been 
obtained in [2J], implying strings may pass through such singularity, the instability of this 
solution was also pointed out in ^j]. Therefore further delicate consideration is needed to 
the understanding of cosmological singularity in string theory. 

Recently another novel program called as discrete gravity has shed new light on this 
problem^] as well. The idea is to discretize Einstein equations. But in contrast to usual 
lattice quantum gravity with fixed multipliers which breaks diffeomorphism invariance mani- 
festly, in discrete gravity the lagrangian multipliers have to vary with "time" level n such that 
discretized equations can be preserved consistent. Furthermore, unlike ordinary mechanism 
of canonical quantum gravity, discrete gravity contains an attractive feature of constraint 
free such that many conceptual problems like the problem of time might be solved in a more 
promising way[(|. In context of cosmology this program has also appealing advantages of 
evolving through the Big-Bang singularity in a general fashion and naturally introducing a 
notion of time at quantum mechanical level. 

In this paper we apply the discrete techniques to study the cosmological model which 

describes the universe may evolve from a big crunch to a big bang. In this context we will 

see our results support the argument that the cosmological singularity can be avoided in a 

rn 

general fashion. But with a distinct starting point from 5], we prefer to considering the 
fate of singularities in universes with A = because in the case of A ^ 0, the classical 
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singularities in vacuum are often coordinate singularities in the absence of matter, although 
this claim does not mean that the singularity at quantum mechanical level must not genuine. 
The discrete theory of FRW cosmology with A = has been considered inj^]. However, as 
its authors pointed out, the discrete theory describing the universe with A = often suffers 
from the absence of the continuum limit, which has negative implications to the quantization 
of the theory. We will also meet this trouble in next section when the model with massless 
scalar field is considered. One purpose of our paper is to see whether we may improve the 
continuum behavior of the discrete theory with zero cosmo logical constant. As a way coming 
out, we consider the universe dominated by Chaplygin gas rather than the ordinary massless 
scalar field in section four. Such an exotic fluid behaves as matter dust at early time, but 
contribute a nonzero effective cosmological constant at later time. 

We organize the paper as follows. In next section we briefly review the framework of 
discrete mechanics focusing on constrained system. Then we turn to consider the model 
with a massless scalar field with A = in section three. After discussing the relations 
between solutions in discretized theory and those in continuum field theory, we argue that 
the singularity may be avoided in discrete cosmology. In section four we study Chaplygin 
gas dominated universe to obtain a better continuum limit of such universes with vanishing 
cosmological constant at early time. 



II. DISCRETE THEORY OF CONSTRAINED SYSTEM 

In this section we briefly review the classical framework of discrete theory of constrained 
system which we will use in next sections. For more details and in particular the quantization 
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We start with the action of a continuum mechanical 



of discrete system we refer to 
system 

S = J dtL c (q,q). (1) 

The discretized formalism is implemented by splitting time into equal intervals t n = ne 
where e is an infinitesimal constant. As a result, all the derivatives with respect to time are 
replaced by discrete first order finite difference, 

q = q n , q = . (2) 
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The time integral in the action then is replaced by a sum S = ^2^=1 L(qn, Qn+i), 
where L(q n ,q n+ i) is understood as the discretized Lagrangian at instant n with relations 
L(q n ,q n+1 ) = eL c (q,q). 

For a constrained system, we can straightforwardly obtain L(q n , q n +i) from L c (q, q) with 
the form, 

L(q n ,q n+1 ) = L(n,n + 1) 

= p n (q n +i - q n ) ~ eH(q n ,p n ) - X n C(q n ,p n ), (3) 

where H is the discrete version of Hamiltonian and C denotes the constraint. Note in the 
last term the parameter e has been absorbed into A n . In this form we have also defined the 
canonically conjugate momenta as 

dL(q n ,q n+1 ) 

Pn+i = — g- • (4) 

In terms of canonical conjugate pairs the equations of motion can be written as, 

_ dH(q n ,P r t +1 ) _ dC(q n ,PZ +1 ) 



dq n " dq n 

dH(q n ,P« +1 ) dC(q n ,P« +1 ) 



qn+i — q n — e TTFfi h A ™ apq 

Ur n+l OJr ri+l 

C(q n ,PZ +1 ) = 0. (5) 

One significant difference between these equations and those continuum ones is that in 
discrete theory the lagrangian multipliers have to be fixed at each instant in order for the 
discrete version of the "constraints" to be preserved. 

Therefore given a continuum theory with constraints we may obtain its discrete La- 
grangian and equations of motion. Going further along this approach, we have two impor- 
tant issues to consider. One is to investigate whether the discrete theory may give rise to 
any different prediction comparing with the continuum theory. The other one is to check 
whether the continuum theory can be recovered by taking the continuum limit of the dis- 
crete theory. In next sections we will apply this framework to cosmology and see these 
cosmological models may help us to clarify these questions. 
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III. BIG CRUNCH/BANG UNIVERSE 



Consider a general 4d theory of a scalar field coupled to gravity. The Lagrangian has the 
form 

S = [ (RV=9- V^0-V(0)), (6) 

where we use the ( — ,+,+,+) signature for spacetime. In context of cosmology we adopt 
minisuperspace description. The metric reduces in conformal coordinate system to 

ds 2 = a 2 (t) [-N 2 (t)dt 2 + dx 2 ] , (7) 

where a is the scale factor and N lapse function. We strict our attention in four dimensional 
spacetime with vanishing V(4>) in this section such that the lagrangian has the form 

C = ^ r (4d 2 - 1 -a 2 <p 2 ). (8) 

The cosmological solution from big crunch to big bang has the form [lj], 

a(t) = a{0)\t\i (j){t) = (f)(0) + Vsin\t\, (9) 

where t is from — oo to +oo. Obviously this solution has a singularity at a = and <fi = — oo 
as t goes to zero. Such a singularity implies that the lagrangian breaks down near the 
neighborhood of that point. As in the quantum or semiclassical theory of hydrogen whose 
apparently continued spectrum is actually discrete at microscopic level, we propose the simi 
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lar consideration should be applicable to the semiclassical theory of spacetime geometry. 
As a result we conjecture a discretized version of lagrangian (JSJ) probably is a good starting 
point to better understand the microscopic origin of the universe. The classical solution 
© may be only valid macroscopically, but should be achieved by taking continuum limit of 
solutions in the discretized theory. 

First consider the canonical formalism of the theory. We define conjugate momenta to 
scalar factor a(t) and (f)(t), 

dC -8a dC 2a 2 

na = ^a- = —' ^ = ^J = W (10) 

while the conjugate momentum of lapse function N, tc^, is zero. The hamiltonian can be 



constructed as 



H = n a a + - £ = -^(^ - ^f) = NH. (11) 

4 4 a 2 
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The discrete version of an action is obtained from © 



L(n, n + 1) = 7r0 n (0 n+ i - <j> n ) + ir an (a n+1 - a n ) 
4 1 4 a 



+ (12) 



where e has been absorbed into the multiplier N n . 

From (j3J) one find the momenta at instants n and n + 1 respectively are 

P n +1 = ^4>n Pn = 7I "0« 

— U -fn — -l<Pn+l - <Pn) H ^2 

pa pa _ 3AU '" 

- 1 n+1 "an n "an 2 fl 3 

N 7T 

r n+l — U r n _ — l tt n+l ~~ a n) g 

^ +1 = = 1(^-^1). (13) 

+ + tt n 

Taking (a n , n , P£ + x, P^+x) as the fundamental variables yields the following compact equa- 
tions 



p<t> p<t> _ 

r n+X r n ~ 









pa pa 

r n+X r n 


2«n 




3A^ n 7T^ n 


■ - 4>n = 


2a n 




1 "n " an 


fln+1 — a n = 


8 


(^n + l) 2 _ 


3(^ + i) 2 
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(14) 

We see is a conserved quantity independent of level n. For convenience we may set it as 
a positive constant C. Plugging it into the last equation, we have 

P„V = (15) 

Then from the other equations in (|14|) it's straightforward to obtain the recurrence relation 
of the scale factor as 



0"n ®n—X 



q, (16) 
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where q can be any positive constant. Given a fixed q and initial conditions, we have two 
kinds of solutions corresponding to the different signs of P^+i hi (fT5j) . 

a n ~ a q n , (fi n ~ O + 2v/3n(g - 1) 

4gggKg-l) /-C 



or 



y + 2V3n(l - gO 



One remarkable result from these solutions is that the lapse N n is determined at each instant 
by the equation, contrasting to the continuum theory where N is a free Lagrangian multiplier 
and may be set to one uniformly Unfortunately we notice from (jlfij) that unless the factor q 
is very close to one the continuum limit can not be achieved as the "time" level n approaches 
to infinity. In particular the continuum solution can only be achieved with a certain accuracy 
as we choose the desired q. Nevertheless, this discrete solution may give rise to the classical 
solution given in 0. To see this we need set the lapse function to unity as n approaches to 
infinity, otherwise the comparison is not meaningful. Consider the solutions (fTTjl and note 
that 

eN ~ N n ~ 4an{a ^~ <ln \ (19) 
we take the continuum limit and set N — 1, leading to 

aa ~ 1, (20) 

whose solution is a ~ t 1 ^ 2 . Here we have also defined that t = ne and a = a " +1 ~° w . If 

e 

q ~ e 6 where 5 is a small number, then we have 5 ~ ^p^, yielding the classical solution 
4> ~ 0o + \/3ln(t). Of course in this sense q is not a strict constant any more but approaches 
to one as n goes to infinity. 

The cosmological solution from big crunch to big bang can be obtained by combining the 
solution (fTHj) with q' < 1 when n < and (fTTjl with q > 1 when n > 0. To be single valued 
at n = 0, we have adjoint conditions ao = a ,C = C',(po = (p' and g + = 2. But we 
note that the conjugate momentum of a n has a jump from 7Tq_ = 2y / 3^ to /TqV = — 2v / 3^-. 
Physically it means the universe contracts as the time level n increases from — oo until n = 



7 



and then rebound at the minimum size of the universe, turning into an expanding phase as 
n > 0. In this discretized version, we find the singularity as t approaches to zero is always 
avoided unless we set ao = by hand. But in generic points olq should be determined by 
the evolution of scale factor from n < to n = 0, which in classical limit corresponds to the 
evolution from t = — oo to t = 0. As a matter of fact, all these intuitive considerations might 
be invalid as the gravitation effect becomes so strong that the quantum theory of gravity 
should be taken into account. Then the quantum states of the universe is the superpositions 
of such states describing expanding and contracting universes. Nevertheless based on this 
heuristic analysis we believe the evolution of universe becomes more delicate due to the 
quantum effect of gravity and rather than a genuine crunch occurs, it might rebound once 
a minimum size of the scale factor is saturated. 

We realize once the quantization of discrete theory is concerned, we need a probability 
explanation of such special initial conditions, namely why the recurrence constant q is close 
to one but not else. In classical theory we are short of such strategy to present a satisfy- 
ing explanation. Based on above consideration, we see the closest argument appearing in 
cosmological model from big crunch to big bang is the requirement that q + q' = 2. If a n 
is always positive as it should be, then q and q' could be any positive number in (0,2), not 
directly leading to our expected value e s and e~ s respectively. One might expect that the 
classical solution (fTBj) is not stable unless the q approaches to unity, in other words q = 1 
might be thought of as an attractor and such running of q from any initial constant to unity 
might be implemented or required at quantum mechanical level. But following perturbative 
analysis shows it's not the case. Consider the perturbation 

a n = a$ + 5a n , (j) n = <f>W + 5(p n , 

p: = p^ + 5p:, pt = p*® + 5 pt 

N n = N®+6N n . (21) 
Substituting into equations of motion (jl4j) and keeping the first order correction, we have 

5a n+1 H — 5a n + q 2 5a n „ 1 = 0. (22) 

On the other hand, perturbations directly based on (fTBj) leads to 

Sa n+1 - 2q5a n + q 2 5a n _ 1 = 0. (23) 



The consistency requires that q ~ 1, which looks satisfying. But another kind of reasonable 
perturbation is a n = ao(q + 5q) n such that Sa n+ i = B ^-q5a n . It immediately leads to the 
conclusion that solutions (fTTj) and (fT%|) are classically stable except for q — 1, at that point 
Sa n could be any value, which is opposite to what we desire. 

We hope this dilemma may be solved by further progress in the quantum theory of gravity. 
But before that we want to present another model where a good continuum limit can be 
achieved nearby the singularity in next section. 



IV. CHAPLYGIN GAS 



Chaplygin gas as an exotic fluid has been greatly investigated in FRW cosmological 
modelp, 3 ( also see |a] for more references). Since it behaves like a cosmological constant 
at later time of the universe, namely contributing a positive energy density and negative 
pressure, it is supposed to be responsible for the acceleration of the universe, as indicated 
by recent observation 0. But at early time of the universe it behaves just like a dust matter 
such that the effective cosmological constant is zero. One interesting feature of Chaplygin 
gas is that its lagrangian is quite similar to the effective lagrangian of tachyon field in string 
theory, which has been proposed by A. Sen[10], even though the exact expression is not clear 
yet. The lagrangian of both Chaplygin gas and tachyon fields can be given as 




£ = ^K + — 3— f-isvl- (24) 

For Chaplygin gas, the potential V(T) is a constant, while for tachyon fields V(T) has a 
maximum a ; r l a,d exponent d™ to^geTQ. 

In this section we consider the discrete theory of Chaplygin cosmology. To compare our 
result with those in previous section, we still adopt the conformal coordinate system 

ds 2 = a 2 (t)[-N 2 (t)dt 2 + dx 2 }. (25) 

As a matter of fact, we are free to redefine the lapse function, 

N'(t) := N(t)a(t), (26) 

and then turn to 

ds 2 = -N'\t)dt 2 + a 2 (t)dx 2 . (27) 



The equations of motion and constraints have the following familiar form 



1 d 
'aWdi \N' 



ItiG V(T) 



1 - 



T 2 

W 1 



3T 2 
2N> 2 



a 2 



8irG V{T) 



a 2 N' 2 
N' 



d It 



AT'2 _ T 2 dt \ N' 



f2 

ht 1 dV n 
+ 3 ^ + VdT=°- 



(28) 



As is familiar, we call the first two equations as Raychaudhuri and Friedman equations 
respectively in standard cosmology and the third one is the equation of motion for tachyon 
fields. Furthermore, as these three equations are not independent we are free to set the lapse 
function N' — 1. 

Next we consider the canonical formalism of the theory. Define the conjugate momenta 
of scale factor a and scalar field T 

dC —a 



7T a 



71t 



da N 

dC 4ttG a 3 VT 



8t 3 y/NW-t 2 ' 



(29) 



then the Hamiltonian is given as 



H = TT a d + 7T T T-C 



-1 / . 2 _ 8nGa 4 N 2 
2N \ a 3 



V 



1 - 



T 2 

W^ 1 



(30) 



The discretized version of the lagrangian is 



L(n, n + 1) = n an (a n+1 - a n ) + ir Tn (T n+1 - T, 



N n 



7T„ 



(31) 



The canonically conjugate momenta of these variables at instant n + 1 and n are respectively 



10 



given by 



pa 
r n+l 



7To 



r n+l — U 

P nT — 



P T 



p N - n 



where we have defined 



A T 



P a = tt 4- TV 

J n — "an T iv n 



47rGV 2 3a° 



P n° = - [K+l - + ^an] 

N n a n 7r Tn 



pTTT 

n 



Pi 



(T n +i — T n ) — 



Xr 



T^Tn 



P N = — (tt 

n 2 



/ AnGal y N n a n V dV 
V 3 J X n dT 

— 2a n X n ), 



V 3 



(32) 



(33) 



In discrete theory the time evolution is described by a canonical transformation from the 
variables from (q n ,p n ) to (q n ,p n+ i). Taking (a n , T n , P£, as the set of basic variables, we 
rewrite these equations of motion as 



pa _ pa _ _nr 



AnGV 2 a e n 



T n +i — T n 



P T - P 1 



N n a n TT T n 



Xr. 



AnGalV N n <LV dV 



A„ dT 



[P: +1 ) 2 - 2a n X n = 



(34) 



Above consideration is applicable to both Chaplygin gas and tachyon field. From now 
on we only focus on the case of Chaplygin gas which has a constant V. The immediate 
consequence from above equation is that P^+i is conserved, which is nothing but 7r Tn . First 
we consider the expanding universe with very small scale factor comparing with Ti Tn , 

ATxGaiV., 



-), 



(35) 
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then A n ~ txtu = Pn+i an d the equations of motion reduce to 

a n+l ~ a n — ~N n 7T an 
P^l - = -NnTTTn 

7n+i — T n = N n a n 

pT _ pT _ q 
r n+l r n — u 

(P n a +1 ) 2 - 2a n P n T +1 = 0. (36) 

After solving for N n and P^ from the second and fifth equation, we may obtain the 
recurrence relations of a n from the first equation as 

a 2 n+ i - Qa n+ ia n + 9a 2 n - 4a„a„_i = 0. (37) 

As we mentioned the evolution scheme is given by the canonical transformation from 
n to n + 1. Given the initial conditions, the value of scale factor can be determined at 
any instant by the recurrence relation (|37|). In particular, whether the universe expands or 
contracts depends only on these initial data. For instance if we take a± > do, then we will 
see a n+ i > a n for any n > 1, implying an expanding universe. On the other side, if we 
choose a -i > a Q , then we find a n+1 < a n for any n < — 1, hence implying an contracting 
universe. Similar to the discussion in previous section, the universe from big crunch to big 
bang may be obtained by connecting these two solutions at n = 0. We may also plot graphs 
to illustrate the discrete evolution of a by choosing appropriate parameters such that the 
scale factor goes to zero as n approaches to from both n < and n > 0, implying a 
singularity at t = in the continuum limit. But in discrete theory more detailed analysis 
indicates that scale factor takes a small but non-zero value at n = in general. Thus its 
solutions may also be used to describe a universe from big crunch to big bang but avoiding 
the singularity. 

We remark that it is a conjecture that the universe from big crunch to big bang may be 
obtained by connecting two solutions at singularity, but not a simple solution to classical 
equations In our discrete theory it says that the set of (a_i,a ,ai) does not satisfy the 
classical equations or ([37)1 . The dynamical explanation of how the universe rebounds 
at the minimum size ao is absent in this classical framework, but it is not unreasonable as 
we expect the quantum theory of discrete gravity would play an essential role to solve this 
issue. 



12 



Now we turn to consider the continuum limit as n goes to infinity. This can be done by 
recovering the first and second order derivatives with respect to time through 



€ 

ttn+i ~ 2a n + a n _i 



(38) 



Then equation (|3*7)l becomes 



a 2 - Aad = 0. (39) 
The general solution to this equation is 

a = c 1 {t + c 2 )K (40) 

where c\ and C2 are two constants, viewed as a remnant of the reparameterization invariance 
of the continuum theory Now going back to the other equations in (J36J) and taking the 
continuum limit we straightforwardly obtain the asymptotical behavior of other parameters 

vr an ~t 2 / 3 N n ~T 1/3 T~f 2 , (41) 
and thus the metric has the form 

ds 2 ~t*/\-r 2 '*dt 2 + dx1). (42) 

To see Chaplygin gas may be though of as a dust matter during this epoch we turn to the 
following conventional coordinate system with lapse function N' — 1: 

ds 2 ~-dT 2 + T*/ 3 dx 2 . (43) 

Therefore the "effective" scale factor behaves as a' ~ T 2 / 3 , indicating an expanding but 
decelerating universe as T > 0. One interesting point here is that during the origin of the 
universe the Chaplygin gas field may be understood as the proper time of the universe. Such 
kind of behavior is something like the time reversal of tachyon field, as pointed out in ^l] , 
where tachyon may play a role of time variable as rolling down to infinity at later time of 
the universe. 

To see our discrete scheme actually leads to the same physics as the continuum field 
theory of Chaplygin gas, we may also consider the continuum limit for large scale factor. In 
this case we have, 

VTTn < ( g )• (44) 
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Since V is a constant, then 



and the equations reduce to 



A„ = (45) 



a n+ i - a n = -N n ix a 
-AN n . 

N n a n 7l T n 



P: +1 - P: = -4iV n A n 



An 

pT _ pT _ q 
r n+l r n — u 

(F n a +1 ) 2 - 2a n F n r +1 = 0. (46) 

Similarly we obtain the recurrence relation of a n 

2a n a n+ i — 3a n + a^_ x = 0. (47) 

Taking the continuum limit yields the differential equation 

2ad + a 2 = 0. (48) 

One solution to this equation is a ~ t 2//3 . Substituting it back to (|3BJ) we obtain the 
asymptotical behavior of other variables as follows 

vr an ~t 4/3 N n ~r*t* r~r 2 . (49) 

As a result, the metric becomes 

ds i^--L dT 2 + T -2/3 dx l (50) 

We note that the Chaplygin gas does not play a role of proper time any more. In effect in 
the coordinate system with N' = 1 the metric goes to 

ds 2 ~ — dr 2 + exp -rdx 2 . (51) 

It indicates the universe is accelerating, thus with a positive cosmological constant, which 
is consistent with the result of continuum theory of Chaplygin gas. 
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V. CONCLUSIONS 



We have applied discrete gravity theory to investigate the universes from big crunch to 
big bang, mainly concerning the following two issues. One is on the fate of cosmological 
singularity. We proposed an attractive point of view on this issue based on the consistent 
discretization approach to general relativity. On the other hand we investigated the contin- 
uum limit of discrete theories describing the universe with vanishing cosmological constant. 
It's clear one may take the attempts presented here as the first step to understand what 
role the discrete gravity would play in quantum gravity and cosmology. In particular in 
this paper we leave all the considerations on the quantization of discrete gravity for further 
investigation. We expect the approach to discrete general relativity will show more insight 
into the problem of the cosmological singularity and quantum gravity. 
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